THE HEIGHT FLUCTUATIONS OF AN OFF-CRITICAL DIMER MODEL 

ON THE SQUARE GRID 



SUNIL CHHITA 

Abstract. The dimer model on a planar bipartite graph can be viewed as a random surface 
measure. We study these fluctuations for a dimer model on the square grid with two different 
classes of weights and provide a condition for their equivalence. In the thermodynamic limit 
and scaling window, these height fluctuations are shown to be non-Gaussian. 



1. Introduction 

In recent years, two dimensional grid based statistical mechanical models have been the 
focus of much research. These discrete models often possess a critical point, where there is 
an abrupt change of behavior, called the phase transition. At this point, the models can 
exhibit macroscopic randomness. This leads to studying the scaling limit, i.e. the behavior 
of the model when the lattice spacing tends to zero. 

Before taking the scaling limit, consider the off-critical model which is a perturbation 
from the critical point. These types of models, in general, do not feature any macroscopic 
randomness and taking the scaling limit of such a model gives a trivial measure. However, 
one could consider the behavior when the lattice size tends to zero while simultaneously 
moving the parameters towards the critical point. We call this limit a scaling window. 
By choosing an appropriate re-scaling, we often find non-trivial scaling window measures 
for off-critical statistical mechanical models. Studying these measures sometimes provides 
additional information about the critical point behavior. 

In this paper, we focus on the dimer model (for the original paper, see [6]) on a bipartite 
planar graph, which is a particular two dimensional statistical mechanical model, defined as 
follows for a finite graph. Call a dimer configuration to be a subset of edges in which each 
vertex is incident to one edge. The dimer model is a probability measure on the set of all 
dimer configurations of the underlying graph. Weights can be applied to the edges so that 
the probability of a configuration is proportional to the product of the edge weights. 

1.1. Overview of Setup. Let G be an infinite planar bipartite graph. As G is bipartite, 
we can define a height function which is in one-to-one correspondence (up to height level) 
with each dimer configuration (see [3). This gives the viewpoint that the dimer model is 
the graph of a random function from — )• M, i.e. a random surface (see [5l [21] for more 
details of random surfaces) . A precise definition of the height function for the dimer model 
can be found in [TT]. Informally, the height function can be defined as follows. 

Suppose d is the number of edges incident to each vertex and that d is fixed for G. The 
height function is a continuous function from to M which is integer valued at the center 
of each face of G. The height change between the centers of two faces sharing an edge is 
ziz{d — 1) if there is a dimer covering the shared edge or =pl if there is no dimer covering 
that shared edge. We choose the following sign convention: suppose that the black vertex 
is to the left when traversing the shared edge between two faces. The height change along 
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the traversal is {d — 1) if there is a dimer covering the shared edge and -1 otherwise. With 
this setup, it can easily be seen that the total height change at the centers of the faces 
surrounding each vertex is zero. 

Define the toroidal graph G„ to be the quotient of G/nZ^, that is, the action of nZ^ on G 
by horizontal and vertical translations. Let A4{Gn) denote the set of all dimer coverings on 
Gn- Define a probability measure Hn for M £ A^(G„) so that the probability of observing M 
is proportional to the product of all the edge weights. For dimer models on periodic graphs, 
there is a limiting measure ^ when the system size is sent to infinity. This limiting measure 
is a Gibbs measure and is called the thermodynamic limit. 

As mentioned above, the scaling window represents the limiting measure (of an observable) 
when the lattice spacing tends to zero while simultaneously moving the parameters to the 
critical point. We use the following notation for scaling windows: for an order e perturbation 
from criticality, the scaling window (e",e^) of a measure is the limiting measure when re- 
scaling the horizontal axis by e", the vertical axis by and e — t- for a,/3 > 0. In other 
words, a box consisting of l/e^"^^) lattice points is sent to a box with continuum area 1 under 
the scaling window (e",e^). In the literature, scaling windows have appeared in different 
guises such as critical windows or intermediate phases. In Section 11.31 we briefly survey 
some of these results for other two dimensional discrete models and provide some motivating 
remarks. 

1.2. Results. This paper considers two diff'erently weighted dimer models of the square grid. 
We say that the dimer model is flipped if the edge weights alternate between ri and r^ along 
a vertical line and between r2 and r4 along a horizontal line. Let /^/(ri, r2, ra, r4) denote 
the flipped dimer model with edge weights ri , r2 , r^ and r4 . We say that the dimer model is 
drifted if the edge weights are si, . . . S4 for edges surrounding vertices Zg, the even sub-lattice 
of Z^. Let fid{si, S2, S3, S4) denote the drifted dimer model with edge weights si, S2, S3, S4. 

Theorem 1.1. The measures /i/(ri, r2, r3, r4) and Hdisi, S2, S3, S4) are equivalent in the 
thermodynamic limit when si = ri/r4, S2 = r|/(rir4), S3 = r|/(rir4) and S4 = r^/vi. 

Using the KPW-Temperley correspondence [13], the drifted dimer model is in bijec- 
tion with random directed spanning trees with the root at infinity. The random span- 
ning tree can be constructed using Wilson's algorithm [21] with random walk probabilities 
(si, S2, S3, S4)/(si + S2 + S3 -|- S4) for the north, east, south and west steps. The flipped 
dimer model also defines a spanning tree (and dual tree) with the root at infinity but this 
tree cannot be constructed using Wilson's algorithm with a massive random walk. It follows 
from Theorem 11.11 that these two trees are equal in distribution. It is not known whether 
there is combinatorial construction which gives a bijection between these two trees. 

The dimer model /i/(l, 1,1,1) (or //(^(l, 1, 1, 1)) is a critical statistical model [8\. We 
consider a perturbation of this dimer model with flipped weights and study the behavior 
of the height fluctuations in the scaling window. As a consequence, this provides results 
regarding the height fluctuations of the drifted dimer model (in the scaling window). We 
find that 

Theorem 1.2. In the thermodynamic limit, the height fluctuations for the flipped dimer 
model fif (1,1 — Aie, 1 — A2e, 1) in the scaling window (e, e) are non-Gaussian. 

The proof of non-Gaussian fluctuations is based on showing that the fourth moments of 
the height function cannot be written as the sum of a product of paired moments, i.e. does 
not satisfy Wick's Theorem. Note that in the limit as Ai and A2 — >■ 0, the fluctuations are 
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Gaussian and are given by the Gaussian Free Field [lOj. If Ai and A2 — > 00, the fluctuations 
are given by White Noise. The probabihty distribution of the fluctuations of the height 
function of /^/(l, 1 — Ai, 1 — A2e, 1) under the scahng window (e, e) remains an open problem. 

I. 3. Related Results and Motivating Remarks. We provide some of the (recent) results 
that are closely connected with scaling limits and windows of the dimer model and random 
directed spanning trees. 

Richard Kenyon, in \W] and [HI, showed that the scaling limit of the height fluctuations 
of the critical dimer model on the square grid is conformally invariant and is given by a 
Gaussian Free Field (see [22j for more details on the GFF). This is the same measure if we 
take the scaling window (e, e) followed by the limit Ai, A2 — )• for ///(1, 1 — Aie, 1 — A2e, 1). 
Recently, there have been many different results such as [20] proving that the contours of the 
discrete Gaussian free field converge to SLE(4) (Schramm Loewner Evolution) in the scaling 
limit. SLE(4) is also the conjectured distribution of paths formed from the scaling limit of 
the critical double dimer model. Excellent surveys of SLE can be be found in [15] and [23] . 

The uniform spanning tree has been of considerable interest and many results can be 
found in [2j. The KPW-Temperely bijection [13] shows that the uniform spanning tree is 
in bijection with the critical dimer model. This bijection allows some further computations, 
which are accessible using the dimer model technology. Indeed, [9] provided the first rigorous 
computational evidence that the uniform spanning tree has some conformal structure. It 
has been shown that the scaling limit of the Peano curve of the 2D uniform spanning tree is 
SLE(8) m)- 

Scaling windows of other discrete models have also been considered. In the literature, 
these have often been called near-critical or intermediate phase. For 2D site percolation on 
a triangular lattice, [18] classified the different scaling windows while [1] shows that there 
are an infinite number of small loops for near-critical percolation under a certain choice of 
scaling window. 

There exists scaling windows of models whose scaling limits give SLE. These objects are 
called Massive SLE [17]. They can be constructed by perturbing the fugacity and performing 
the scaling window with respect to the fugacity. For example, the off-critical LERW features 
loops with a small weight and the lattice is then re-scaled while simultaneously sending the 
small weight to zero (see [1]). It is not known whether the scaling windows of the dimer 
models studied in this paper are connected with Massive SLE. 

Non-Gaussian fiuctuations are not uncommon for 2D statistical mechanical models per- 
turbed from criticality. For example, the so-called 2D (/)^-model has non-Gaussian fluctua- 
tions in the scaling limit [3]. 

The motivation behind Theorem 11.21 was to flnd properties of the height fluctuations 
between the scaling limit of the critical and off-critical dimer model on the square grid. 
In order to determine whether the off-critical model under re-scaling provided some non- 
trivial behavior, we needed to send the parameters of the model to the critical point while 
simultaneously sending the mesh size to zero which is deflned as the scaling window in Section 

II. 11 The choice of the scaling window parameter exponents a and /3 are chosen so that the 
correlation length of height function is finite under re-scaling. Other symmetric choices of 
the scaling window parameter exponents leads to the fiuctuations governed by either the 
Gaussian Free field or White Noise depending on whether the re-scaling is too slow or too 
fast. 

Finally, many physical models are studied away from the critical point using perturbation 
theory with the goal that they belong to some universality class. For an example related to 
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the dimer model on the square grid, [13] considers the six- vertex model close to the critical 
point and shows the scaling under a certain limit is the same as the critical sine-Gordon 
model [16] . We believe that the model studied in this paper belongs to the universality class 
of free field theory deformed by a certain operator but the exact details of this remains open. 

1.4. Overview of the Paper. The paper is organized as follows. In Section [21 we provide 
some background information about the dimer model, such as, the Kasteleyn method. The 
proof of Theorem 11.11 the bijections of the discrete models as well as the operators for the 
underlying random walks are given in Section [3l In Section U we compute the relevant 
Green's functions for the operators, which allow us to prove Theorem II. 2i In the same 
section, we also provide the amoeba interpretation of the scaling window. 

2. Setup 

2.1. Partition Function, Local Probabilities and Notation. We refer to the graph Gi 
as the fundamental domain. Orient the edges of Gi so that the number of counter clockwise 
edges per face is odd. This is called the Kasteleyn orientation. Assume that G (the infinite 
planar graph) has the periodic Kasteleyn orientation induced by Gi. As G is periodic, we 
can write any vertex of G as a vertex of Gi plus a translation of the fundamental domain, 
where the first co-ordinate of the translation is in the direction of the vector (1,0) and the 
second co-ordinate is in the direction of the vector (0,1) (i.e. v + {x,y) is the vertex v in 
the fundamental domain translated by {x,y)). For example, the fundamental domain of 
the square octagon lattice with the diagonal edges assigned weight t and off diagonal edges 
assigned weight 1 is depicted in Figured] 




Figure 1. Fundamental Domain of the square octagon lattice. Take Kaste- 
leyn orientation with the convention that all the edges are directed from white 
vertices to black vertices except for the edges (^1,62) and (1^3,64) which are 
both oriented in the reverse direction. 

Let B and W represent the white and black vertices of a bipartite graph. For any bipartite 
graph, define a weighted adjacency matrix by K : WxB — )• R, with K{w, b) given by the weight 
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and direction of the edge wb. This matrix is called the Kasteleyn Matrix. For translational 
invariant graphs, it is easier to work with a Fourier transform as this gives information about 
the free energy and provides a tractable way of computing edge probabilities. We now setup 
the Kasteleyn matrix on the torus using the method outlined in |12j . 

Let V{Gi) = {wi,bj} be the set of vertices of Gi labelled as per Figured) Let 71 (and 
resp. 72) be a path in the dual of Gi, winding around the torus in the direction of the 
vector (1,0) ((0,1)). Write eij for the weight of the edge Wibj. For parameters z and w, let 



eij, where 



(2.1) 



ni(n2) 



/ 1 




Wjb^ crosses 71(72) from below 
does not cross 71(72) 
y — 1 Wibj crosses 71(72) from above 



Define Ki to be the weighted adjacency matrix of Gi with entry given by the weight 
connecting the i*'' white vertex to the j^^ black vertex and sign determined by the Kasteleyn 
orientation. This means. 



(2.2) 



{Ki{z,w)) 







if Wi ~ bj , and Wi 
if Wi ~ bj , and bi 
otherwise 



Wi 



where Wi ~ bj means that Wi and bj share an edge and Vi — )• Vj denotes an arrow from Vi to 
Vj in the Kasteleyn orientation. For the above square octagon example, we have 



(2.3) 



Ki{z,w) 
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t 
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t 
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w z 
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The relevance of the above discussion is that we can write the partition function and 
probabilities of observing a local configuration in the thermodynamic limit as expressions 
involving Ki{z,w) (see [12j). 

The characteristic polynomial for the fundamental domain, P{z,w), is defined to be the 
determinant of Ki{z,w) (see [12j). For the square octagon example, this is given by 



(2.4) 



P{z,w) = 4 + 



w 



We can write the free energy as the following: 



(2.5) 



logZ 



(47r2) 



dw dz 
[ogP[z, w) 

--lJ\w\ = l W Z 



We can now give the local statistics formula [7J which is given by the inverse Kasteleyn 
entries. Let E = {ei, . . . , e;} be a set of edges with Cj = {w*, b*). 



[Wj,b*, 



det ({K-\blw*: 



(2.6) Fieu...,em)= \^^K{t 

where, assuming b and w are in the same fundamental domain. 



(2.7) 



K \b,w + ix,y)) 



(27ri) 



[Z,W)hwW 



/l<i,j<m 



,dw dz 



W Z 
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2.2. Gauge Transformation and P-equivalence. We say that two weight functions, vi 
and V2, are gauge equivalent if there are functions on the white vertices and black vertices, 
/i and /2, such that for each edge e = wb, vi{e) = fi{b)f2{w)v2{e). The mapping from vi 
and V2 is cahed a gauge transformation and is measure preserving. 

We say that two dimer models with periodic weights are P-equivalent if the edge weights 
of the local configurations of each dimer model differ by a multiplicative constant dependent 
on the locations of each local configuration. It is not necessarily true that P-equivalence 
preserves the measure. 

An example of P-equivalence is the following: consider a dimer model on the square grid 
which has edge weights for edges incident to the vertices (2x, 2y), edge weight for the 
edge connecting the vertex (2x — 2,2y-\- 1) to the vertex (2x — l,2y + 1) and edge weights 
for the rest of the edges incident to the vertices {2x — l,2y + 1) for x,y € Z and s ^ 1. This 
dimer model is P-equivalent to the dimer model on the square grid with each edge having 
weight s as shown in Figure [2j Note that this is not a gauge transformation. 



-Q 



-Q 



-Q 



Figure 2. An example of P-equivalence. 



2.3. Amoebas and Phase Diagrams. The amoeba of P{z,w), denoted A(P), is defined 
to be the set {(log |z|,log \w\) : P{z,w) = 0} which forms curves in M^. These curves split 

into different regions, with each region describing the phase behavior of the measure of 
the height function for the dimer model, fi, with characteristic polynomial P{z, w), on some 
bounded domain. For the height differences can be deterministic (frozen), stochastic with 
unbounded variance (liquid) or stochastic with bounded variance (gaseous). It has been 
shown [12], that n is frozen, liquid or gaseous, in the closure of the unbounded component 
of M2\A(P), in the interior of A(P) or in the closure of the bounded component of A(P). 
The origin provides the specific location of the measure. A more in depth discussion can be 
found in [121 H]. 

An example of an amoeba the dimer model on the square octagon lattice, ()2.4p . can be 
found in Figure O For t < ^/2, the phase diagram features solid, liquid and gaseous phases 
and the measure is in the gaseous phase. At t = \/2, there is no gaseous phase and the 
measure is in the liquid phase. The height fluctuations in the scaling limit are given by 
a Gaussian Free Field ( [HI EJ [l2j ) . For t = \/2 — e, the height fluctuations in the scaling 
window (e, e) represent the transition between the liquid and gaseous phases. 

2.4. Some Notation. For the dimer models on the square grid, we have the following 
notation: Let W and B denote the white and black vertices. Let Bq denote the black vertices 
{(23;, 2y) : x, ?/ G Z}, Bi denote the black vertices {(2x -|- 1, 2y -|- 1) : x, y E Z},Wo denote the 
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Figure 3. The amoeba of the dimer model on the square octagon lattice 
with characteristic polynomial (I2.4p . with phases labelled. 



white vertices {{2x + 1, 2y) : x,y £ Z} and Wi denote the white vertices {(2x, 2y + 1) : x,y £ 
Z}. 

3. The Equivalence of /xj and fi^ 

In this section, we prove Theorem 11.11 By setting si = ri/r^, S2 = r2/{rir4), S3 = 
r'^/irir^) and S4 = r^/ri, the measures r2, rs, r^) and fJ-disi, S2, S3, S4) are P-equivalent, 

as shown in Figure [H 



-O- 



-O 
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S3 
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1 



Figure 4. P-equivalence between the flipped dimer model to the drifted 
dimer model with si = ri/r^, S2 = r2/{riri), S3 = r'^/{riri) and S4 = r^/ri. 



Using the Fourier Transform method, the characteristic polynomials for both models can 
be found and are given by 

(3.1) w) = r2r4 [w + w~^) + rir^ [z + z'-^) - {rf + r2 + r| + r|) 
for the flipped dimer model and 

(3.2) Pd{z, w) = S2W + S4tt;~"^ + S32; + siz"-^ - (si + S2 + S3 + S4) 




Figure 5. The amoebas for the flipped dimer model (left) and the drifted 
dimer model (right). The horizontal and vertical lines are the x and y axis. 

Proof of Theorem \l.l[ Let bi E Bj and Wi G Wj for i G {0,1}. The fundamental domain 
for the square grid and the Kasteleyn orientation are given in Figure [6j Let Kf(z,w) and 
Kii(z, w) be the Kasteleyn matrices for the flipped and drifted dimer models, where the rows 
are indexed by wq and wi and the columns are indexed by bo and hi. These are given by 

(3.4) Kj{z,w) = ( ^^-J^r, 

^ \ — r\w —r2 + r^z 

and 

(3.5) Kf{z,w) = _i -1 , _i 

\ S3 — SlW —l + z 

where the signs are determined by the Kasteleyn orientation and the appropriate edge weights 
can be found in Figured] 
We have that 

(3.6) Kd{z, w) = Di.Kf (^z, "^w) .D2 



where 

(3.7) Di -. 
and 

(3.8) D2 



1 

riT4 





r2 ri 


r-3 



rs 




rir4 



Let a = r4/r2 and /3 = r^/ri and let Kj:^ and ^ represent the flipped and drifted 
inverse Kasteleyn matrices respectively. Suppose that Wi{x,y) represents Wi in the {x,y) 
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Figure 6. The Kasteleyn orientation and vertex labels for the square grid. 

fundamental domain from hj for i, j S {0, 1}. In order to distinguish between the vertices, 
write Kj^{bi,Wj{x,y)) for the inverse Kasteleyn entry between bi and Wj{x,y). We have 

(3.9) K^\bo,wo{x,y)) 

(3.10) K^\bo,wiix,y)) 

(3.11) K^Hbi,Mx,y)) 
and 

(3.12) K-\b^,wi{x,y)) 

It remains to check that any configuration of edges is equivalent using the local statistics 
formula (j2.6p . The above values have the property that for i,j G {0, 1}, 

(3.13) Wkie)Kj\bi,Wjix,y)) = Wd{e)a'' K^\h,Wj{x,y)) 

where e is the edge incident to bi in the same fundamental domain, Wk and Wd are the edge 
weights of e in the flipped and drifted dimer models. 

For the plane, each entry in the inverse Kasteleyn matrix for the flipped dimer model is 
the same as the drifted dimer model up to a multiplicative factor dependent on the distance 
of the two vertices (j3.13p . Each combination of edges is given by a determinant of the inverse 
Kasteleyn matrix. However, the sum of the exponents (in the drifted dimer model) for each 
term in the expansion of the determinant is zero. Hence, the determinant expansions for the 
flipped and drifted dimer models are the same. Therefore, for any A €z J-, fJ-f^A) = fj,d{A). 
This means that on the plane fJ-f = almost surely. 

□ 

An immediate corollary of Theorem 11.11 is 

Corollary 1. The height function defined from fif has the same distribution as the height 
function defined from Hd ■ 



n 

a^'p^+y 
n 

a-+^f3y . 



Kj'^{bo,wo{x,y)), 
Ky'^{bo,wi{x,y)), 
Ky^{bi,wo{x,y)), 



{bi,wi{x,y)). 
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3.1. Bijections with other Models. Denote ST(n, e, s, to be the (random) directed 
spanning tree with weights n = n{yi,y2) upwards, e = 6(7/1,1/2) eastwards, etc. These 
weights can be dependent on the position of the edge. The language of spanning trees is 
similar to the language of the dimer model and so to avoid ambiguity, we refer to the vertices 
and directed edges of the directed spanning trees as nodes and arrows. 

The drifted dimer model, /^d(si, S2, ss, •S4) is equal in distribution with ST(si, S2, S3, S4). 
The flipped dimer model, /^/(n, r2, ra, r4), is equal in distribution to another spanning tree 
model, which is gauge equivalent to ST(s^^^^^, Sj^^^^, sl^"*"^^, 54^^^^) which is P-equivalent 
to ST(si,S2, 53,84). 

Let nsoit) denote the measure of the dimer model on the square octagon lattice (with a 
specified boundary conditions) with edge weights t on the diagonal edges and edge weight 
1 on the remaining edges. The case t = \/2 corresponds to the dimer model on the square 
octagon lattice being at criticality. Let t = 

Lemma 3.1. /i5o('S^/2) with any boundary conditions is equal in distribution to /i/(l, s, s, 1) 
and P-equivalent to ;Urf(l, s^, s^, 1) and ST(1, 1, s^, s^), with the appropriate boundary condi- 
tions. 

Proof. Most of the proof follows from the corresponding result for the uniform measure in 
|13j . However, we do provide details in obtaining the appropriate weights of the spanning 
tree. The first step is to decompose the square octagon lattice to a square grid using urban 
renewal. For the square octagon lattice labelled in Figuredl allocate the vertices bi,wi, 62, 1^2 
to be rich cities and squares with the vertices 63, ws, 64, to be poor cities. We can decom- 
pose each poor city along with the connecting edge and vertex of each of the four neighboring 
rich cities into a new city with a different weight, say s, along each edge in the new city as 
shown in Figure [71 




Figure 7. The urban renewal transform, which transforms a poor city 
(neighbored by 4 rich cities), into a new city, with weights s on the edges 

For the edges surrounding a poor city, there are either 0, 2 or 4 of the edges covered by 
dimers. In the decomposition by urban renewal, these correspond to no edges, 1 edge or two 
edges covered by dimers in the new city. As we have set t = s^/2, it is clear that the weights 
of the local configurations agree. 

We color in blue the upper left hand vertices of each new city. We color in red the lower 
right hand vertices of each new city. The red vertices form the set of nodes for the spanning 
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tree while the blue vertices form the nodes of the associated dual spanning tree. In order 
to obtain the appropriate weights on the arrows of the spanning tree, we must first apply 
a gauge transformation so that all the edges incident to a blue vertex have weight 1. An 
instance of this gauge transformation is shown in Figure [8j The gauge transformed dimer 
model has different weights depending on the locations of the edges, that is, the dimer 
model is ;L(rf(s-(j'i+2^2)^ 5-(yi+j/2+2)^ ^-(4/1+1/2+2)^ ^-(yi+j/a)), it ig dear that this dimer model is 

P-equivalent //^(l, s^, s^, 1). 



1 1 



1- 1 



1 



Figure 8. The gauge transformation, which applies a transformation on the 
edge weights. Here, v = s^. 



-O- 



-O 



□ 

The significance of the colored vertices in the above proof is that the blue vertices are the 
vertex set of the directed spanning tree and the red vertices are the vertex set of the dual 
spanning tree. A dimer covering a blue vertex and a white vertex corresponds to a directed 
edge from that blue vertex passing through the white vertex to the next blue vertex. 

The urban renewal transformation guarantees that the height function defined on the 
dual of the rich cites is the same as the height function on the dual of the square grid. More 
precisely. 

Lemma 3.2. The urban renewal transformation is a measure preserving bijection for the 
height function. 

Proof. The height function is defined up to height difference which is determined by the 
presence of a dimer in the configuration. As the urban renewal transformation is measure 
preserving, then it remains to check that the height change across the edges are the same. 
This can be validated by checking each local configuration. □ 

3.2. Massive and Drifted Operators. In this subsection we introduce the massive and 
drifted discrete operators. Applying these operators to the inverse Kasteleyn entries gives a 
discrete partial differential equation with the boundary conditions given by the graph of the 
underlying dimer model. For f G W, define the Kasteleyn operator by 

(3.14) Kfiv) = J2K{v,u)f{u) 

where n ~ u means that u and v are nearest neighbors, and (f , uY^ entry of the K{u, v) is 
the Kasteleyn matrix. For v £ B, define the dual Kasteleyn operator by 

(3.15) K*f{v) = Y,Kiv,u)fiu) 
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Let ei = (2,0) and 62 = (0,2). The massive operator is given by 
Lff = K*Kf{v) 
(3.16) 

= r2rif{v + ei) + r2rif{v - 62) + rir3f{v + 62) + rirzf{v - 62) - (r^ + r| + r| + r|). 

Define the shifted graph of 1? to be the graph on J? with the weights of the horizontal 
edges shifted upwards by one lattice spacing and the weights of the vertical edges shifted to 
the right by one lattice spacing, as seen in Figure [H For ti G B, define the dual Kasteleyn 







1 
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^1 



Figure 9. The original graph and the shifted graph for the flipped dimer model. 



operator on the shifted graph by 

(3.17) kf{v) = Y,K{v,u)f{u) 

where K{v, u) is the (f , uf^ entry of the Kasteleyn matrix on the the shifted graph. The 
drifted operator is given by 

Ldf{v) = KKf{v) 

(3.18) = sif{v + 62) + S2f{v + ei) + S3f{v - 62) + s^fiv - ei) - (si + S2 + S3 + S4). 

4. Scaling Window 

This section focuses on the behavior of the height function from the flipped dimer model. 
In doing so, we have to introduce the Green's function of the massive random walk as well 
as the Green's function of the drifted random walk. These functions allow the computation 
of the inverse Kasteleyn entries in the scaling window, which can be used to show that the 
fluctuations in the scaling window form a non-Gaussian fleld. 

4.1. Green's Functions. In this subsection, we compute the massive and drifted Green's 
functions in the scaling window. These functions are closely related to the flipped and drifted 
dimer model ~ the inverse Kasteleyn entries for the flipped (and respectively drifted) dimer 
model can be written as a linear combination of massive (drifted) Green's functions. 

4.1.1. Drifted Green's Function. The drifted operator, defined in ()3.18p . can be re-written 
as 

(4.1) Ldf{v) = sidyf{v) + S2dxf{v) - s^dyfiv - 62) - s^d^fiv - ei) 

where dx and dy represent the discrete derivatives. Let be the Green's function of the 
on the plane, where the superscript represents the lattice spacing. In other words, we have 

(4.2) -LlH}=6o. 
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There is a finite drift in the scaling window (e,e) if si — S3 and S2 — S4 are both 0(e). 
Choosing si and S4 equal to 1, we have 

Lemma 4.1. Suppose that linie-s-o ^/^i^s ~ ^i) = A2 > and iim^^o l/e(s4 — S2) = Ai > 0. 
In the scaling window (e, e), for {x,y) € M^\{(0,0)} 

(4.3) H^aix, y) = ^ exp (Aix - Asy) Bk (o, ^(Af + A2)(x2 + y2)^ 

where Bx{n, z) is the modified Bessel function of the second kind with index n. 

Proof. For the above choice of weights, the operator is the Markov generator of a two- 
dimensional random walk with drift (Aie, — A2e) because 

(4.4) L'^Jiv) = A^fiv) + A2e(/(^ - ei) - f{v)) + Xieif{v - 62) - f{v)) 

where A*^ is the discrete Laplacian operator on eZ^. By Donsker's theorem, this random 
walk with drift converges (weakly) to a Brownian motion with drift A := (Ai, — A2) and the 
partial differential equation now reads 

(4.5) {-A + X-V)H'a=5o. 

As a two-dimensional Brownian motion with drift is transient, the Greens' function can be 
computed using the usual probabilistic methods, namely, for y E M^\{0} 

(4.6) H^aiO,y) = / PiO,y,t)dt 



(4.7) =/ ^e.J-\y^]dt 



2TTt ^ V 2t 

where p{0, y, t) is the transition kernel for the Brownian motion with drift. The integral in 
(j4.7p is a standard integral and is exactly the Bessel Function given in the lemma. 

□ 

The function is the Green's function for a two dimensional Brownian motion with drift 
on the plane. This is the reason behind the name - 'the drifted dimer model'. 

4.1.2. Green's function of the Massive Random Walk. The operator defined in (j3.16p is the 
(Markov) generator of a massive random walk. The random walk has transition probabilities 
(rir3,r2r4,rir3,r2r4)/(r^ + '^2 + '^i + '^4) the north, east, south and west steps and 
probability ((ri — r-^)'^ + (r2 — r4)^)/(r^ + '^i + '^s + ^4) visiting the so-called graveyard. 
The massive operator (j3.16p can be rewritten as 

(4.8) Ljf= {r2ridl + r^r^dl - (ri - r^f - (r2 - r^f) f 

where dx and dy represent the discrete derivatives. Let Hj be the Green's function of on 
the plane, where the superscripts represents the lattice spacing. In other words, 

(4.9) - L}H} = 60. 
Note that we can also write 

w^z^ dz dw 



(4.10) H}{x,y) = -^[ [ 



,|=iy|2|=i Pf{z,w) z w 

i.e. Hj is the Fourier Coefficients of 1/ Pf{z,w). This can be seen by evaluating the appro- 
priate linear combination of Lj. 
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In order for the massive random walk to have a finite rate of visiting the graveyard in 
the scaling window (e, e), we require ri — and r2 — r4 to be ce where c is some constant. 
Choosing ri and r4 equal to 1, we have 

Lemma 4.2. Suppose that lime_!.o l/e(?'3 — ri) = Ai > and lime_>.o l/e(?'4 — T2) = A2 > 0. 
In the scaling window (e, e), for {x,y) G M^\{(0,0)} 

(4.11) Hjix, y) = ^Bk (0, ^(A2 + Ai)(x2 + y2)^ 

where BK{n, z) is the modified Bessel function of the second kind with index n. 

Proof. Set Gi{m,n) = H^{m,n)r'^r^'^ for {m,n) £ 1? . G\ saftisfies the discrete partial 
differential equation ()4.9p and so G\ = Hj . is computed in Lemma 14.11 and evaluating 
the exponential factors gives Hj. 

□ 

If ri and r4 do not approach 1 in the scaling window but ri — r^ and r2 — r^ are 0(e), then 
the Green's function in the scaling window does exist. It is given by a stretch or contraction 
in the appropriate direction of the Green's function given in Lemma 14.21 

4.2. Non-Gaussian Scaling Window. In this subsection, we show that the height func- 
tion for the flipped dimer model with specific edge weights converges to a non-Gaussian 
field, i.e. we prove Theorem 11.21 In order to prove this assertion, we need to determine the 
correlation functions for the flipped dimer model in the scaling window. Theorem 11.21 follows 
by showing that the fourth moment of the height function cannot be written as the sum of 
a product of the second moments of the height function and hence violates Wick's theorem. 

The next lemma gives the explicit expressions for the correlation function of the dimer 
model. 

Lemma 4.3. Let |A| = a/A| + A^, z = (x,y) and \z\ = x^ + y^, with Uiq^^q l/e{bi — Wj) = 
{x,y), where hi G Bj and wj G Wj for i,j £ {0, 1}. In the scaling window (e, e), the entries of 
the inverse Kasteleyn matrix are given by 



(4.12) 


lim —K 

e^O e 




wo) 


= %xBk (1, |A||z|) + -XiBk (0, |A||z|) , 


(4.13) 


lim —K~ 




Wl) 


= -^^yBK (1, |A||z|) + -X2BK (0, |A||z|) 
Tr\z\ vr 


(4.14) 


lim —K 

e-S>0 e 




,Wo) 


= %yBK (1, \m) + -^2Bk (0, \X\\z\) 

TT\z\ TT 


and 








= ^xBk (1, |A||z|) - -XiBk (0, |A||z|) . 
TT\z\ vr 


(4.15) 


lim —K~ 

e-S>0 e 




■Wl) 



where Bxin, r) is the modified Bessel functions of the second kind with index n. 

Proof. The entries of the inverse Kasteleyn matrix are given by the local statistics formula 
(j2.6p where 

(4.16) k{z,w) = ( '^""1 ''"'-1 V 

^ ' ^ ' y r^ — z —r2 + w J 
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Each entry of the inverse Kasteleyn matrix can be written as a hnear combination of the 
massive Green's function (in the discrete setting). For example, with lattice spacing e, 
K~'^{bo,wo) is given by 

(4.17) -r2Hl{x,y)+H'a{x-e,y). 

Using a Taylor expansion and Lemma [4.2^ this can be rewritten in terms of Bessel functions. 
For example, K^^{bQ,wo) now reads 

(4.18) ,(^.^H^^^x,y) + X,HUx,y)^ +0(6^). 
Dividing through by e and taking limits gives the results. 

□ 

Let Sn{vi, . . . , Vn) denote the n point correlation function for the height change along the 
edges fi, . . . , u„ in the scaling window (e, e). Suppose that 71, . . . , 7^ are n disjoint lines in 
the scaling window (e, e), we let 

(4.19) 5„(7i, . . . ,7„) = / ... / Sn{vi,...,Vn)dvi...dvn- 

J-yi J "in 

Let 71, . . . , 74 be four lines (in the scaling window) on the same vertical line with 7in7j = 0. 
Any two point correlation function can be determined for the lines 71 , . . . , 74 

Lemma 4.4. In the scaling window (e, e), the two point correlation function for the height 
function along 7^ and 7j for i ^ j and i,j £ {1, 2, 3, 4} is given by 

(4.20) 52(7i,7j) = ^IAP / / BKiO,\XH-Vj\f + BKil,\X\\vi-V2\fdv^dvj. 
Proof. The definition of the height function gives 



(4.21) 52(7i,7i) = lim32E 



2n 2m 
fc=0 l=Q 

where the index k and I are the discretized index of the lines 7^ and 7^-, and is the 
indicator event that there is a horizontal dimer crossing 7j at index k. Interchanging the 
sums, we are required to compute 

(4.22) E [{X2i - X2i+i){X2, - X2j+i)\ . 

This quantity is a linear combination of four 2 by 2 determinants, i.e. 

(4,23, ti-^'''^^{ K-C,M^ ) 

n,ii=o 

where lime_j.o l/e|'U^m — bn\ = \m — n\. These entries can approximate by Lemma 14.31 This 
gives 

I \ 12 

(4.24) - 2e^H^ {Bk{0, \X\^\v, - v,\f + Bk{1, \Xf\vj - Vi\ f) 

where vi = k + 0{e) and V2 = / + 0(e). Plugging (I4.24p back into (I4.21j) and taking the limit 
as e — 7- 0, the Riemann sum converges to the required integral. 

□ 
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In a similar fashion, we can find an exact expression for the four point correlation function 
for the change of heights along 71, . . . ,74. However, this expression is rather complicated, 
and the full expression is not required for the proof of Theorem II .21 Instead, we just require 
the following 

Lemma 4.5. In the scaling window (e, e), the four point correlation function for the height 
function along 71 , . . . , 74 is given by 

4 

(4.25) 5*4(71,..., 74) - ^ J|5'2(7i,7j) = / ... / f(x^M)i^'i-'X2,X3,X4)dx4^. . .dxi. 

where /(Ai,A2)(^1' 2:3, 2:4) is defined in the Appendix equation IIA.I]) . When (Ai,A2) > 0, 
we have f(x-,,x2)ixi,X2,X3,Xi) / 0. When (Ai, A2) 0, we have f(^Xi,X2)i^^^ ^"^^ ^3, Xi) 0. 

Proof. From the definition of the height function, we can evaluate 5'4(7i, 72, 73, 74), using 
the four point correlation function of the dimer model, namely. 



(4.26) lim 3^E 



(4.28) det 



2n 2m 2p 2q 

''^{^^t — X2i+l){X2j — X2j+l){X2k — X2k+l){X2l — -'^2/+l) 
i=0 j=0 k=0 1=0 

where the index i,j,k,l are the discretized index of the lines 71 , . . . , 74 , and Xr is the indicator 
event that there is a horizontal dimer crossing the corresponding vertical line. We can 
approximate 

(4.27) E [{X2i - X2i+l)(X2j - X2j+l){X2k - X2k+l){X2l - X21+1)] 

using Lemma 14.31 and then follow the same approach given in Lemma 14.41 i-S- we can write 
1 

^_iyi+ji+ki+h 

ii,ji,ki,h=0 

( K-\wi,,bj,) K-\w,,,bk,) K-\wi„bi^) \ 

K-^{wj^,bi^) K-^{wj^,bk^) K-^{wj^,bi^) 

K-^{wk^,bi^) K-^{wki,bj^) K'^{wk^,bi^) 

\ K-\wi,A,) K-\wi„bj,) K'\wi,,bk,) I 

with lime_j.o l/e|t(;mi ~ ^ml = \'m — n\. Each of the entries of the above matrix can be 
approximated using Lemma 14.31 and each term contains a factor or e. We can expand out 
the above determinants and (I4.26P is a Riemann sum which converges in the limit e — )> to 
a sum of products of Bessel functions. 

In order to prove Lemma 14.51 it remains to subtract the sum of the products of the 
two point correlation function. This leaves the integral of some linear combination of Bessel 
functions which are exactly the contributions of the 4-cycles in the above 4 by 4 determinant, 
i-*^- !^^--- !^^f{\i,\2)i^^^^2,X3,XA)dxi...dxi. The full expression of /(Ai,a2)(2;i, 2:2, X3, 2:4) 
can be found in the Appendix, equation (jA.ip . Each term of /{Ai,A2)(^1' ^2, a^s, 2:4) is of the 
form 

(4.29) ± 4^5x(ii, \\\zi)BK{i2, \\\z2)BK{i3, | A|z3)Si^(i4, |A|z4) 

where each ik G {0, 1} for 1 < < 4, zi + Z2 + ^3 + «4 is divisible by 2 and for 
1 < s < r < 4 . Notice that /(Ai.A2)(^1) ^2, 2:3, X4) is a function of A/Af + Af and is non-zero 
for |A| > 0. 
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Taking the generalized series expansion around |A| = (using the generahzed series ex- 
pansion of the Bessel functions), as Bk{^, z) is approximately l/z + z/A{—l + 2Eg — log4 + 
2\og{z)) + 0{z'^) and Si^(0, z) = log(2;) - + log(2) + 0(2:^) where Eq is Euler's constant, 
we find that the highest order (with respect to | A|) of the integrand is 0(— |Ap log |A|). Using 
the dominate convergence theorem, the left hand side of (j4.25p tends to zero as |A| tends to 
zero. 

□ 



Proof of Theorem Lemma [4. 51 and Lemma [A . 1 1 given in the Appendix shows that we can 
find intervals 71 , • • • , 74 such that 

4 

(4.30) 54(7i,---,74)/ 5]n'^2(7.,7r) 

This means that the correlation functions do not satisfy Wick's theorem (that the moments 
of the normal can be written as the permanent of the Green's function). This proves the 
theorem. □ 

Corollary 2. For \ > 0, in the scaling window (e, e) the height fluctuations of ^so{{X~^^)'^) 
on the plane are non-Gaussian. 

Proof. By Lemma [3. II and Lemma [3.21 the height fluctuations of ^soUX ~ -^^)^) ^^^^ equal in 
distribution to the height fluctuations of /Uj(l, 1 — Ae, 1 — Ae, 1), which are non-Gaussian by 
Theorem O 

□ 

The first part of Theorem 11.21 also holds when ri yA- 1, r4 yA- 1 as e — t- but ri — r^ and 
r2 — r^ is 0(e). If ri — t- k2 and r4 — )■ /ci, the argument in the Green's function becomes 



(4.31) 



The arguments in the Green's function change by a multiplicative factor and hence, the 
derivative terms in the inverse Kasteleyn entries change by multiplicative factors. However, 
the same proof for the first part of Theorem 11.21 holds, but with these extra multiplicative 
factors attached to the derivatives of the Green's function. 

4.3. Amoeba Interpretation. The amoeba for the flipped dimer model /Uj(l, 1 — Aie, 1 — 
A2e, 1), consists of an ellipse with the intercept with the line el*"' = given by 



(4.32) 



iblog 



2 - 2eA2 + (Af + A^) - Je^ (Af + A^) {e^j + (2 - eAs) 2) 



2 - 2eA2 

and the intercept with the line e'^' = given by 

"2 + 2eAi + e2(Af + Al) - Je^ (A? + A^) (^2^2 + (2 + eAi) 2) 



(4.33) 



±log 



2 + 2eAi 



An expansion of these terms gives 



(4.34) ±eJXl + Xl±e^X2\Xl + \l 



18 



SUNIL CHHITA 



and 



(4.35) ±e\l\l + \l±e^\i\\l + \^ 



Rescaling the ellipse (by e) and as e tends to zero, the ellipse tends to a circle with radius 
y^Af + A^. For measures of the form /xj(ri(e), r2(e), r3(e), r4(e)), with lim^^on 7^ linie->.o ?'4, 
the bounded region is an ellipse, with foci 0(e) provided ri — and r2 — r4 is 0(e). 

As the height fluctuations of /i/(l, 1 — Aie, 1 — A2e, 1) and /Urf(l, (1 — Aie)^, (1 — A2e)^, 1) 
are equivalent, then the fluctuations at the boundary of the gaseous region and the center 
of the gaseous region are equivalent in the scaling window. In fact, the height fluctuations 
of the whole gaseous region is equivalent which can be seen by gauge transformations. 
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Appendix A. 



In the appendix, we present the function from Lemma [4.51 For yi = X2 — xi, 2/2 = 2^3 — X2 
and 2/3 = X4 — X3, we have 
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- 4\X\^Bk[0, \X\yi]BK[0, |A|(yi + 2/2)]Sk[0, \X\y3]BK[0, \X\iy2 + ys)] 

- 4\X\^Bk[0, \X\yi]BK[0, \X\y2]BK[0, \X\y3]BK[0, |A|(yi + + ys)] 

- 4|A|^Bk[0, \X\y2]BK[0, |A|(yi + y2)]BK[0, \X\{y2 + y3)]BK[0, |A|(yi +y2 + 2/3)] 
+ 4\X\^Bk[0, \X\y3]BK[0, |A|(yi + 2/2 + y3)]BK[l, \X\yi]BK[l, \X\y2] 

- 4\XfBK[0, \X\ys]BK[0, \X\{y2 +y3)]BK[l, |A|yi]B^[l, |A|(yi + ya)] 

- 4\X\^Bk[0, \X\{y2 + y3)]BK[0, |A|(yi + y2 + y3)]BK[l, \X\y2]BK[l, \X\{yi + ys)] 
+ 4|A|^B,^[0, |A|(yi + y2)]5x[0, |A|(y2 + y3)]BK[l, |A|yi]S^[l, |A|y3] 

- 4|A|^Bk[0, |A|y2]Si^[0, |A|(yi + y2 + y3)]S^[l, \X\yi]BK[l, \X\y3] 
+ 4|A|45k[0, |A|yi]5;^ [0, |A|(yi + y2 + y3)]5K[l, |A|y2]5x[l, |A|y3] 
+ 4\X\^Bk[0, \X\yi]BK[0, |A|(y2 + y3)]Bx[l, |A|(yi + y2)]Bx[l, lAjys] 
+ 4\X\^Bk[0, |A|(yi + y2)]5i^[0, \X\y3]BK[l, \X\yi]BK[l, \X\{y2 + ys)] 

- 4|A|^S,^[0, |A|(yi + y2)]5i^[0, |A|(yi + ya + y3)]5i^[l, \X\y2]BK[l, \X\iy2 + y3)] 
+ 4\X\*Bk[0, |A|yi]5x[0, |A|y3]Si^[l, |A|(yi + y2)]5x[l, |A|(y2 + ys)] 

- 4\X\^Bk[0, \X\y2]BK[0, |A|(yi + ya + y3)]5x[l, |A|(yi + y2)]BK[l, \X\iy2 + ys)] 

- 4\X\^Bk[0, \X\yi]BK[0, |A|(yi + y2)]Sx[l, |A|y3]5x[l, |A|(y2 + ys)] 

- 4\X\^Bk[1, |A|yi]i?x[l, |A|(yi + y2)]i?x[l, |A|ys]i?x[l, |A|(y2 + ys)] 

- 4\X\^Bk[0, \X\y2]BK[0, |A|y3]S,^[l, |A|yi]5x[l, |A|(yi + ys + ys)] 
+ 4|A|4Si^[0, |A|yi]5x[0, |A|ys]S,^[l, |A|y2]Bx[l, |A|(yi + y2 + ys)] 

- 4\X\*Bk[0, |A|(yi + y2)]Si^[0, |A|(y2 + ys)]5K[l, |A|y2]5x[l, |A|(yi + y2 + ys)] 

- 4\X\^Bk[0, \X\y2]BK[0, |A|(y2 + ys)]5x[l, |A|(yi + y2)]5x[l, |A|(yi + y2 + ys)] 

- 4|A|^S,^[0, \X\yi]BK[0, |A|y2]5i^[l, |A|ys]Sx[l, |A|(yi + ys + ys)] 
+ 4|A|^S,^[1, |A|yi]Bx[l, |A|y2]S,^[l, |A|ys]5x[l, |A|(yi + ys + ys)] 

- 4\X\*Bk[0, \X\y2]BK[0, |A|(yi + y2)]Sx[l, |A|(y2 + ys)]Sx[l, |A|(yi + y2 + ys)] 
(A.l) - 4\X\*Bk[1, \X\y2]BK[l, |A|(yi + y2)]5x[l, |A|(y2 + ys)]^x[l, |A|(yi + y2 + ys)] 



We also have the following lemma which is based on a crude argument. 
Lemma A.l. There exists 71, ... ,74, so that 



Proof. To show that J ■ ■ ■ I-y /ai,A2(3J1) 3:^2, iCs, a^4)da;i ... dx4 7^ for some choice of intervals 



7i, choose 71 = (0, 1/2"), 72 = (1 - 1/2", 1), 73 = (2, 2 + 1/2") and 74 = (3, 3 - 1/2"). We 



claim that /ai,A2(^1) ^2, X3, X4) is positive over this interval and so the integral over 71, . . . , 74 



(A.2) 




IS non-zero. 
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We can bound yi and 2/3 in terms of 2/2 where yi,y2 and are defined above. These are 
given by 

^ 2" - 2 

and 

2" — 2 

(A.4) ^2 _^ ^ <yz<y2- 

We can separate (jA.ip into positive terms and negative terms, i.e. we can write (jA.ip as 
2/2) ys) — 5-(yi) 2/2, ys) where g'_ and g+ are positive. A lower bound for the positive 
terms can be achieved by setting yi and 7/3 to 2/2- An upper bound of the absolute value of 
the negative terms can be achieved by setting yi and 2/3 to 1/2(2" — 2)/ (2" + 2). In other 
words, we have 

(A.5) 5+ (2/1, 2/2, 2/3) > 5+ (2/2, 2/2, 2/2) 

and 

, . /2"-2 2" -2 

(A.6) g- (2/1, 2/2, 2/3) <g- ( 2» + 2 ^'^'^^' 2" + 2 ^^ 

The difference between these two bounds is positive provided we take n = |A|m for |A| > 1 
and n = m/|A| for |A| < 1 where m > 20 (these bounds are not tight). 

Indeed, by setting (2" — 2)/(2"' + 2) = 1 — x, we can take a series expansion with respect 
to X between the positive and negative terms, which is given by 

(A.7) 5r+(7/2,?/2,2/2) - y- (2/2, 2/2, 2/2) " a:^!!. (2/2, 2/2 , 2/2) + 0{x^) 

where the derivative is with respect to x. The term is strictly positive and is bounded 
below by e'^^^^IAj^ ^j^jg 

comparing the following positive and negative terms of (jA.ip 
and noting that the difference in each case is bounded below by e~^^^^l'^L These differences 
are given explicitly by 

(A.8) 4|A|4i3i,[0, \\\y2fBK{\ 1\\\y2f - m^^K^ I A|2/2]'i?x[0, 2|A|y2]' 

(A.9) 4|A|4i?K[0, 2|A|y2]'i?A'[l, |A|y2]' - 4|A|4s^[0, \\\y2fBK{^, 3|A|2/2] 

m^BK% \\\y2\BK% 3|A|2/2]i?x[l, |A|y2]' - 4|A|4i?K[0, \\\y2\BK% 2\\\y2fBK% 3|A|y2] 

- %\\tBK% 2|A|2/2]5x[0, 3|A|2/2]5k[1, |A|2/2]Sx[1, 2|A|2/2] 

- m^BK% \\\y2fBK{\, |A|y2]i3i^[l,3|A|y2] 
(A.IO) 

-4|A|^Si^[0,|A|y2]i?K[0,3|A|y2]i?/dl,2|A|y2]' 

4|A|4Si^[l,|A|y2fi3i^[l,3|A|y2]-4|A|"S^[l,|A|y2]'i?i^[l,2|A|y2]' 

-4|A|4i?K[0,2|A|y2]2i?x[l,|A|y2]Si^[l,3|A|y2] 

- 8|A|^fiK[0, |A|y2]fix[0, 2|A|y2]S^[l, 2|A|y2]5x[l, 3|A|y2] 

(A.ll) -4|A|^i?A-[l,|A|y2]i?A'[l,2|A|y2]2i3i^[l,3|A|y2]. 

For |A| > 1, the coefficient of the x^ term in the above series expansion is negative and 
its absolute value is bounded above by e"'^'^^ for y2 > 1. As x is smaller than 1/2^"'^, then 
/ai.A2(^i, ^2, 2^3, 2:4) is positive over 71,..., 74 for |A| > 1. For < |A| < 1, the absolute 
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value of the coefficient of x is bounded above by some constant c > 2 for ?/2 ^ 1- As 
1/2'^^'^ is less than e"^^^^''', the order term is less than the order term. This means 
f\i,X2{xi^X2,X3,X4) > over 71, ... ,74. 

□ 
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